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Individuals with major depressive disorder (MDD) display abnormal neurophysiological responses to 
psychological stress but little is known about their neurophysiological responses to physiological stressors. 
Using [ 15 0]-H 2 0 positron emission tomography we assessed whether the regional cerebral blood flow 
(rCBF) response to arterial cannulation differed between patients with MDD and healthy controls (HCs). 
Fifty-one MDD patients and 62 HCs were scanned following arterial cannulation and 15 MDD patients and 
17 HCs were scanned without arterial cannulation. A region-of-interest analysis showed that a significantly 
increased rCBF of the anterior cingulate cortex and right amygdala was associated with arterial cannulation 
in MDD. A whole brain analysis showed increased rCBF of the right post-central gyrus, left temporopolar 
cortex, and right amygdala during arterial cannulation in MDD patients. The rCBF in the right amygdala 
was significantly correlated with depression severity. Conceivably, the limbic response to invasive physical 
stress is greater in MDD subjects than in HCs. 

I ndividuals with major depressive disorder (MDD) show altered behavioral 1 , neuroendocrine 2 , autonomic 3 , 
I and neurophysiological 4 responses to psychological stress. Neuroimaging paradigms have therefore been 
I developed which implicitly or explicitly attempt to mimic aspects of the social environment in order to better 
model neural responses to psychosocial stress, and thereby improve our understanding of MDD. Stimuli have 
taken the form of emotionally-valenced words or visual stimuli, sad, angry or frightened faces, and recordings of 
critical comments. 

In aggregate, the results of these studies suggest that MDD is characterized by a dysregulation of the viscer- 
omotor network, which regulates endocrine and behavioral responses to stress 5 . In particular, there is evidence for 
a reduction in "top-down" cortical (prefrontal cortex) control over structures involved in mediating the visceral 
and behavioral aspects of emotional behavior, such as the anterior insula (AI), ventral striatum, mediodorsal 
thalamus, hippocampal subiculum, amygdala, hypothalamus and brain-stem nuclei. 

Nevertheless, the functional imaging data are largely based on responses to what are, by practical and ethical 
necessity, relatively subtle and situationally-constrained emotional or psychosocial interventions rather than 
responses to both non-social and physiological stressors outside of the laboratory. One such environmental 
stressor, which has been understudied in patients who meet DSM-IV criteria for MDD, is the response to robust 
physiological stressors involving pain or mildly invasive procedures. 

While conducting PET neuroreceptor mapping studies of MDD, our laboratory routinely acquired ls O-H 2 0 
regional cerebral blood flow (rCBF) scans prior to neuroreceptor imaging to facilitate data processing and/or 
tracer modeling. For some of these studies, radial arterial cannulation is necessitated by the need for an arterial 
input function. This mildly invasive procedure typically induces mild-to-moderate levels of pain (despite the 
application of local anesthetic agents), possibly via stimulation of the radial nerve 6 , and activates the neural 
systems implicated in stress biology, including the limbic-hypothalamic-pituitary-adrenal (HP A) axis 7 , resulting 
in elevated Cortisol 8 , and a vascular sympathetic autonomic response. Since a tight coupling (positive correlation) 
exists between changes in rCBF and local terminal field metabolic activity 910 , the methodological requirements of 
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these neuroreceptor PET studies afforded us an opportunity to study 
the neurophysiological response of MDD patients and healthy con- 
trols to a more invasive and aversive stressor than has been used in 
previous neuroimaging studies. 

The aim of our current study was to test whether neurophysiolo- 
gical regulation of response to an unpleasant physiological stressor is 
modified by MDD. Specifically, we aimed to test whether compared 
to HC subjects, MDD patients show a greater difference in rCBF 
between high and low stress experimental conditions. That is, are 
depressed subjects more sensitive neurophysiologically to an invas- 
ive stressor than healthy controls? 

Stress responses are mediated by the sympathetic-adrenomedul- 
lary (SAM) and the HPA axes, resulting in increased release of 
epinephrine, norepinephrine, and Cortisol, respectively 11 . In two 
independent cohorts we previously found that compared to healthy 
controls, depressed patients show increased Cortisol release in res- 
ponse to arterial cannulation 8-12 . Limbic structures such as the amyg- 
dala, which project directly to the paraventricular nucleus (PVN) of 
the hypothalamus, regulate stress-induced glucocorticoid secretion 
by disinhibiting cortisol-releasing factor (CRF) release from the 
PVN 13 . Conversely, there is an extensive body of research dem- 
onstrating that epinephrine, norepinephrine, and Cortisol influence 
amygdala function 1415 . Based on these data, as well as the extensive 
literature on the role of the amygdala in depression and pain, we 
expected to find increased rCBF of the amygdala in MDDart+ 
patients. 



This study also specifically investigates two regions of interest 
(ROIs) comprising the cortical components of the stress matrix that 
process the emotional, endocrine and autonomic dimensions of 
aversive stimulation - the anterior cingulate (ACC) and anterior 
insular cortices 1617 . In their meta-analysis, Apkarian et al. show that 
45 out of 68 studies of pain induction in healthy subjects reported 
activity changes in insula, the primary viscerosensory and nocicep- 
tive cortex 18 . This meta-analysis also noted changes in ACC activity 
in 46 of 68 studies, consistent with evidence that the ACC modulates 
amygdala- driven CRF release 19 in response to noxious emotional or 
physiological stimuli 18 . 

Results 

Demographic and Clinical Variables. These data are shown in 
Table 1. 

There was no significant difference between the 4 groups in gender 
distribution (x 2 =1.0, p = 0.81) or age at scan (F = 1.9, p = 0.13). No 
significant difference was found between the MDD patients who 
underwent arterial cannulation (MDDart+) and MDD patients 
who did not undergo arterial cannulation (MDDart-) groups on 
the number of weeks medication free (t= 1.2, p=0.24), age at onset 
(t=0.7, p = 0.50), or illness duration (t=1.2, p = 0.25). 

Both the art+ and the art- MDD groups scored higher than the 
art+ and art- HC groups on the MADRS (F=215, p<0.001) but 
there was no significant difference between the art+ and art- MDD 



Table 1 | Demographic and Clinical Characteristics of the Sample 



Trait 




Diagnosis 








Statistics 






MDD_art 


MDD_no_art 


HC_art 


HC_no_art 


Statistic 


DF 


p-value 


Result 


N 


51 


15 


62 


17 










Gender 










x 2 =1.0 


3 


0.81 


NS 


Males 


19 


6 


28 


6 










Females 


32 


9 


34 


1 1 










Percent Female 


63 


60 


55 


65 










Age at Scan 


37.6±9.7 


3 1 .9± 1 0.7 




34.8±9.6 


31.9±8.4 


F= 1 .9 


3 


0.13 


NS 


















Past Medication Use ( 3 weeks 


















previous; N) 






NA 


NA 










Naive 


20 


10 














AD 


27 


5 














Mood Stabilizers 


7 


2 














Benzodiazepines 


2 


1 














Weeks Medication Free 


93.4± 120.7 


27.7±9.3 


NA 


NA 


t=1.2 


32 


0.24 


NS 


Age at Onset 


18.0±6.2 


19.8±1 1.7 


NA 


NA 


t=0.7 


48 


0.50 


NS 


Years ill 


17.8±1 1.4 


13.7±9.7 


NA 


NA 


t=1.2 


48 


0.25 


NS 


Family History of Affective Illness (N) 


















Present 


31 


1 1 


0 


0 










Absent 


4 


1 


62 


17 










Uncertain 


16 


3 


0 


0 










Past History of Substance/ Alcohol 


















Abuse (N) 






NA 


NA 










Present 


8 


2 














Absent 


43 


13 














Common Comorbid Disorders in 


















Sample (N) 






NA 


NA 










Panic Disorder 


5 


3 














PTSD 


4 


4 














GAD/Social Anxiety Disorder 


3 


2 














Eating Disorder 


2 


1 














Mood Ratings 


















MADRS 


24.7±8.1 


22.5±5.4 


0.3±0.7 


0.5±1.2 


t= 1.0° and 


53 a and 


0.34° 


NS 












t=0.7 b 


68 b 


and 0.49 b 





Note: (a) comparison of MDDart+ and MDDart- on the MADRS. 
(b) comparison of HCart+ and HCart- on the MADRS. 
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groups (t=1.0, p=0.34) or between the art+ and art- HC groups 
(t=0.7, p=0.49) on the MADRS. 

Although no demographic or clinical variables differed signifi- 
cantly across the 4 subject groups, the MDD art+ subjects were 
on average 5.5 years older than the MDD art- subjects. Similarly, 
the HC art+ group was on average 3 years older than the HC art- 
group. To ensure that our results were not systematically biased by 
age effects, we therefore added age as a covariate in the general 
linear model. 

Imaging. ROI Analysis. Contrast A: MDDart+ versus MDDart-. 
MDDart+ patients displayed higher rCBF compared with 
MDDart- patients in the following regions (table 2 and fig. 1A): 
the supragenual and pregenual ACC, left posterior insula, and left 
and right amygdalae. 



MDDart+ patients displayed lower rCBF compared with 
MDDart- patients in the following regions (table 2 and fig. IB): left 
and right anterior insulae and the sgACC. 

Contrast B: HCart+ versus HCart-. Within the predefined ROI no 
areas were identified where rCBF differed significantly between the 
HCart + and HCart- groups. 

Contrast C = A-B. The MDD patients showed a significant increase 
in rCBF in response to arterial cannulation in the right amygdala and 
the supragenual ACC relative to healthy controls (table 3, fig. 1C). No 
significant decrease in the rCBF of MDD patients in response to 
arterial cannulation stress was observed relative to healthy controls. 

In post hoc assessments the flow obtained over the locus where the 
peak voxel-t-value was evident within the right amygdala (table 3) 
was significantly correlated (positively) with the MADRS score in the 



Table 2 | Results of th 


e Region 


of-lnterest (ROI) analysis for contrast A showing areas of increased (top panel) or decreased (lower panel) 


rCBF in MDDart+ su 


bjects versus MDDart- 


subjects 








CONTRAST A -ROI Analysis 










Increase 


d rCBF in MDDart+ 


vs. MDDart- 








FDR 








Ouster Size T score 


Z score 


p-value 


p-value 




Coordi nates" 


Region of Interest'' Additional Comments 










X 


Y Z 




2508 4.41 


4.26 


<0.001 


0.006 


-2 


26 24 


ACC, supragenual Cluster size corrected p-value 














<0.001; FWE corrected 














p-value=0.018 


4.25 


4.12 


<0.001 


0.006 


2 


32 15 


ACC, pregenual Cluster size corrected p-value 














=0.002; FWE corrected 














p-value=0.031 


4.00 


3.88 


<0.001 


0.006 


2 


15 25 


ACC, supragenual 


120 3.85 


3.75 


<0.001 


0.006 


24 


1 -20 


R amygdala 


232 3.83 


3.73 


<0.001 


0.006 


-40 


-12 -6 


Insula, L posterior insula 


3.52 


3.44 


<0.001 


0.007 


-44 


-19 5 


Insula, L posterior insula 


72 3.39 


3.32 


<0.001 


0.008 


-24 


1 -20 


Lamygdala 


159 3.39 


3.32 


<0.001 


0.008 


40 


-17 -14 


Insula, R posterior insula 


3.33 


3.26 


0.001 


0.009 


42 


-1 -12 


"Insula", R temporopolar 














cortex 


3.10 


3.04 


0.001 


0.013 


46 


9 -1 1 


"Insula" R temporopolar 














cortex 


CONTRAST A - ROI Analysis 










Decreased rCBF in MDDart+ 


vs. MDDart- 














Additional 


Cluster Size T score 


Z score 


p-value 


FDR p-value 


Coordinates" 


Region of Interest'' Comments 












X Y 


Z 


385 4.51 


4.35 


<0.001 


0.014 




36 -1 


15 "Insula", R dorsal insula/ FWE corrected 














frontal operculum p-value = 0.014, 














Cluster corrected 














p-value = 0.036 


3.53 


3.45 


<0.001 


0.036 




32 22 


10 "Insula", R dorsal insula/ 














intrasulcal lateral OFC (BA47s) 


320 3.41 


3.33 


<0.001 


0.028 




-40 -7 


17 "Insula", L dorsal insula/ Cluster corrected 














frontal operculum p-value = 0.053 


3.17 


3.1 1 


0.001 


0.043 




-40 18 


6 "Insula", L dorsal insula/ 














frontal operculum 


3.04 


2.99 


0.001 


0.050 




-30 8 


1 1 "Insula", L Anterior Insula 


34 3.06 


3.01 


0.001 


0.048 




12 25 


10 Right ACC, subgenual 


a. Coordinates correspond toth 


e stereotaxic 


array ofTalairacha 


id Tournoux 55 


and are in mm from the stereotaxic orig 


n (anterior commissure) with positive x = right, positive y = anterior, and positive z — 


dorsal to a plane containing both the anterior and posterior commissures. 








b. The region-of-interest in whic 


l differences 


n normalized flow 


s listed first; however, because the regional template for the insula extended outside the cortex of primary interest to encompass the entire 


operculum, where the coordinates localized cortical areas other 


than the insula based on 


he stereotaxic coordinates (Mai et al. 54 ; Talairach and Tournoux 55 ; Ongur et al. 20 ) these cortices are named 


accordingly and the ROI name appears in quotation marks (i.e. 


"insula"). 








Significance threshold for peak voxel in cluster: psO.001 . 










Significance threshold for contiguous voxels 


p<0.01 and region > 100 mm 3 


in size. 






Abbreviations: ACC= anterior 


cingulate cortex; L= left; OFC=orbitofrontal cortex; R=right; FDR=false discovery rate corrected p-value; FWE=family-wise error corrected p-value. 
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1 



2 



3 



4 






v=i 




y = -l 



z = 15 





y = 2 



z = -20 



x = 0 



Figure 1 Image sections from a statistical parametric map of voxel t-values which show regions where rCBF is significantly increased or decreased in 
contrasts A and C of the ROI analysis. Stereotaxic coordinates provided beneath each section image denote the distance in mm from the stereotaxic 
origin (anterior commissure), with positive x indicating right of midline (for sagittal sections), positive y indicating anterior (for coronal sections), and 
positive z indicating dorsal (for horizontal sections). Abbreviations: A=anterior to anterior commissure, P=posterior to anterior commissure, L=left 
hemisphere, R= right hemisphere, spACC=supragenual ACC, pgACC=pregenual ACC, AMY= amygdala, INS = insula, sgACC=subgenual ACC. Row A 
shows regions of the brain that show greater rCBF in MDDart+ versus MDDart- patients (contrast A). Al: sagittal slice showing areas where voxel t- 
values show significantly increased activity in the supragenual and pregenual ACC [x=— 2, y=26, z=24). The posterior ACC is also evident but is below 
the significance threshold; A2: coronal slice showing higher activity in the amygdala, bilaterally, and in the supragenual ACC (x=24, y= 1, z= — 20); A3: 
axial slice showing elevated activity of the posterior insula, bilaterally (x= —40, y= — 12, z= — 6). Row B shows regions of the brain that demonstrate 
reduced flow in the MDDart+ group relative to the MDDart- group (contrast A). Bl and B2 display: coronal and axial slices, respectively, showing areas 
of decreased rCBF in the dorsal anterior insula/frontal operculum, bilaterally (x=36, y= — 1, z= 15), while B3 displays an axial slice showing decreased 
rCBF in the sgACC (x= 12, y=25, z= — 10) MDDart+ group relative to the MDDart- group. Row C shows regions of the brain associated with greater 
rCBF in response to arterial cannulation in the MDD + group compared with the MDDart- group relative to differences between the HCart + and HCart - 
groups (contrast C). CI: coronal slice showing areas where flow is increased in the amygdala, bilaterally (x=24, y=2, z= —20), C2: axial slice showing 
areas of increased flow of the amygdala, bilaterally (x=24, y=2, z= —20), C3: sagittal slice showing region of peak rCBF in the supragenual ACC (x=0, 
y=27, z=28). A4, B4, and C4: Key showing correspondence between image color and voxel t-value. 



MDDart+ group (r=0.272, p = 0.045). In the left amygdala the nor- 
malized activity showed a nonsignificant trend towards correlating 
positively with the MADRS score (r=0.231, p = 0.076) (fig. 3). There 
was no significant association between the MADRS scores and nor- 
malized activity of the right amygdala (r=0.036, p = 0.449) and left 
amygdala (r=— 0.33, p = 0.454) in the MDDart- group. There was 
also no significant correlation between the normalized activity of 
the peak voxel in the supragenual ACC cluster and MADRS 
score in either the MDDart+ (r=0.002, p = 0.988) or the MDDart- 
(r= -0.254, p = 0.361) groups. 

Contrast D: MDDart- versus HCart-. Within the predefined ROI no 
areas were identified where the MDDart- group showed greater rCBF 



than the HCart- group. The MDDart- group did however show 
reduced rCBF in the left insula compared with the HCart- group 
(table 4). 

Whole Brain Analysis. Contrast A: MDDart+ versus MDDart-. 
Regions in which MDDart + patients showed increased rCBF 
compared with MDDart- participants included (table 5 and 
fig. 2A): the supragenual and pregenual ACC, right entorhinal 
cortex, temporopolar cortex bilaterally, rostral temporal sulcal 
cortex (junction of superior and middle temporal gyri) bilaterally, 
ventromedial PFC, left ventrolateral PFC, left lateral orbitofrontal 
PFC, postcentral gyri and adjacent inferior parietal cortex bila- 
terally, and precuneus. 
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Table 3 | Region-of-interest analysis showing increased rCBF in the MDD sample versus the HC sample in the high versus low stress 
condition (contrast C) a 

CONTRAST C - ROI Analysis 



Cluster Size T score Z score p-value FDRp-value Stereotaxic Coordinates Region of Interest Additional Comments 

X Y Z 

88 3.30 3.23 0.001 0.330 24 2 -20 R amygdala L amygdala significant 

at peak voxel threshold 
of p<0.002 (cluster 
size 5 1 ) 

140 3.17 3.11 0.001 0.330 0 27 28 ACC, supragenual ACC 

Significance threshold for peak voxel in cluster: p<0.001 . 

Significance threshold for contiguous voxels: p<0.01 and region > 100 mm 3 in size. 
Abbreviations and coordinates interpreted as in table 2. 

a. No regions showed statistically significant decreases in rCBF in the MDD sample versus the HC sample. 



MDD art+ patients showed reduced rCBF compared with the 
MDD art- patients in the following regions (table 5 and fig. 2B): right 
anterior insula, left insula, and left frontal operculum. 

Contrast B: HCart+ versus HCart-. Relative to the HCart- group the 
HCart+ group showed higher rCBF in the right premotor cortex and 
lower rCBF in the left parahippocampal gyrus (table 6). 

Contrast C = A-B. The MDD patients showed a significant increase 
in rCBF compared with HC subjects in response to arterial cannula- 
tion in the following regions (table 7, fig. 2C): the right postcentral 
gyrus, right amygdala, and left temporopolar cortex. 

No regions of the brain showed a decrease in rCBF in MDD 
patients compared with HC subjects in response to the arterial line 
stressor. 

In post hoc assessments there was no significant correlation 
between the normalized activity at the peak voxel in the left tempor- 
opolar cortex and the MADRS scores in either the MDDart+ 
(r=0.108, p = 0.508) or the MDDart- (r= -0.084, p=0.766) groups. 
In addition, there was no significant correlation between the normal- 
ized activity at the peak voxel in the right post-central gyrus and the 
MADRS scores in either the MDDart+ (r=0.154, p= 0.344) or the 
MDDart- (r=— 0.208, p = 0.457) groups. The corresponding rela- 
tionship at the peak voxel in the right amygdala was not tested 
because of its proximity to the peak coordinate in the right amygdala 
in the ROI analysis (see below). 

Contrast D: MDDart- versus HCart-. The MDDart- group showed a 
significant increase in rCBF compared with HCart- subjects in the 
right posterior hippocampus (table 8). 

No regions of the brain showed a significant decrease in rCBF in 
MDDart- patients compared with HCart- subjects. 

Discussion 

The main focus of this study was to determine if there is a pattern of 
rCBF that distinguishes MDD patients from healthy individuals 



when they are subjected to an invasive stressor. In general, depressed 
subjects showed neurophysiological responses to arterial cannula- 
tion stress that were exaggerated relative to those of healthy controls 
in limbic/paralimbic regions of the anterior cingulate cortex, amyg- 
dalae and temporopolar cortex. Our most salient results were as 
follows: 

(1) In the supragenual, pregenual ACC, and the ventromedial 
prefrontal cortex, rCBF increased in MDD patients in the 
arterial cannulation versus the control condition (contrast A: 
table 2, fig. 1A; table 4, fig. 2A), while the corresponding con- 
trast in healthy controls (contrast B) revealed no significant 
difference in these areas. These regions form part of a medial 
prefrontal or "visceromotor" network that plays major roles in 
modulating the visceral (autonomic, endocrine) responses and 
processing associated with stress and emotion 20 . Additionally, 
after controlling for differences in rCBF in healthy subjects 
across high and low stress conditions, a more circumscribed 
region of the ACC showed greater rCBF in MDD patients 
compared with HC subjects (contrast C: table 3, fig. 1C). We 
have labeled this region as the "supragenual ACC" to signify 
that it is immediately superior to the genu of the corpus cal- 
losum but a variety of different terms - ACC, rostral ACC, 
anterior midcingulate cortex, perigenual ACC, and pregenual 
ACC- have been used to describe this region in the literature. 
The supragenual/pregenual ACC has been hypothesized to be 
the human homologue of the rat prelimbic cortex which exerts 
a predominantly inhibitory effect on HPA function via its 
connections with the bed nucleus of the stria terminalis and 
the paraventricular nucleus of the hypothalamus 19 , and is 
highly sensitive to psychological 21 and physiological stress 22 . 
In humans, this region is also activated by noxious stimuli, 
particularly visceral as opposed to cutaneous stimuli, as well 
as psychological stressors, particularly those associated with 
the emotion of fear 23 . Critchley et al 24 showed that neuro- 
physiological activity in the supragenual and pregenual ACC 



Table 4 | Region-of- interest analysis for contrast D showing increased rCBF in the MDD sample versus the HC sample in the low stress 
condition (i.e. neither group underwent arterial cannulation)" 

CONTRAST D - ROI Analysis 

Cluster Size T score Z score p-value FDRp-value Stereotaxic Coordinates Region of Interest 

X Y Z 

33 3.35 3.27 0.001 0.420 -43 -26 4 "Insula", L temporopolar cortex 

Significance threshold for peak voxel in cluster: p<0.001 . 

Significance threshold for contiguous voxels: p<0.01 and region 5 100 mm 3 in size. 
Coordinates are interpreted as in table 2. 

a. No regions showed statistically significant decreases in rCBF in the MDD sample versus the HC sample. 
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Table 5 Results of the whole brain analysis for contrast A showing regions of increased (top panel) and decreased (bottom panel) rCBF in 


MDDart+ subjects versus MDDart 


subjects 


















CONTRAST A - 


Whole Brain Analysis 














Increased rCBF in 


MDDart + vs. MDDart 








Cluster Size" 


T score 


Z score 


i 

p-value 


FDR p-value 


Stereotaxic Coordinates 


Region 












X 


Y 


z 




1 1 75 


4.41 


4.26 


<0.001 


0.040 


-2 


26 


24 


Supragenual ACC 




4.25 


4. 1 2 


<0.001 


0.040 


2 


32 


15 


Pregenual ACC 




a r\r\ 
4.UU 


TOO 
J.OO 


<0.001 


0.040 


-2 


13 


31 


Supragenual ACC 


679 


A O O 

4.38 


A 1 A 

A.Z4 


<0.001 


0.040 


22 


7 


-21 


R entorhinal C 




4.03 


3.92 


<0.001 


0.040 


42 


14 


-23 


R temporopolar C 




3.79 


3.69 


<0.001 


0.040 


44 


1 


-17 


R rostral STG/MTG 


551 


4.28 


4. 1 4 


<0.001 


0.040 


-34 


5 


-22 


L temporopolar C 




o.yo 


T QO 
O.OZ 


<0.001 


0.040 


-40 


-5 


-13 


L rostral STG/MTG 




3.73 


3.69 


<0.001 


0.040 


-34 


12 


-34 


L temporopolar C 


61 


4. 1 0 


3.98 


<0.001 


0.040 


-57 


-33 


33 


L inferior parietal lobule 


72 


3.84 


3./2 


<0.001 


0.040 


59 


-25 


34 


R postcentral gyrus 


40 


3.74 


3.65 


<0.001 


0.040 


50 


-13 


-33 


R inferior temporal gyrus 


65 


3.73 


3.64 


<0.001 


0.040 


0 


26 


-25 


Ventromedial prefrontal C 




3.41 


3.33 


<0.001 


0.047 


-2 


40 


-26 


Ventromedial prefrontal C 


22 


3.62 


3.54 


<0.001 


0.042 


-46 


22 


19 


L ventrolateral prefrontal C 


36 


3.57 


3.51 


<0.001 


0.045 


-42 


-12 


39 


L precentral gyrus 


28 


3.52 


3.44 


<0.001 


0.045 


-42 


41 
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L ventrolateral prefrontal C 




3.22 


3.16 


0.001 


0.052 


-40 


45 


-2 


L lateral orbitofrontal C 


58 


3.40 


3.33 


<0.001 


0.047 


2 


-60 


34 


Precuneus 


31 


3.47 


3.40 


<0.001 


0.046 


-20 


38 


26 


L dorsolateral prefrontal C 
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Decreased rCBF in MDDart+ vs. MDDart- 
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T score 


Z score 
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171 


4.51 


4.35 


<0.001 0.177 


36 


— 1 15 "Insula", 


R dorsa 


insula/frontal o 


Derculum 


49 


4.01 


3.90 


<0.001 0.195 


26 


20 6 "Insula", 


R dorsa 


insula/intrasulcal lateral OFC (BA47s) 


22 


3.41 


3.33 


<0.001 0.403 


-40 


-7 17 "Insula", 


L dorsal 


insula/frontal operculum 


a. Note that the difference in cluster 


sizes shown here vis a vis those in table 2 


eflects the constraint of the anatomical extent o 


the cluster by the ROI boundaries 


Only the local maxima with the three highest 


voxel t-values within a given cluster 


are listed here 


so these loci should be considered additive to those listed in Tables 2 and 3 in the manuscript. 


Significance thre 


hold for peak voxel in cluster: p<0.001 , and cluster size c 


f contiguous voxels for which p<0.001 corre 


iponding to > 1 00 mm 3 . 




Coordinates are 


nterpreted as in table 2. 














Abbreviations: ACC - anterior cini 


ulate cortex; C 


- cortex; STG/MTG - cortex along the su 


oerior temporal sulcus at the j 


jnction of su 


perior temporal gyrus and middle temporal gyrus; PFC - prefrontal 


cortex. 



















correlated with heart rate variability across a variety of mildly 
stressful cognitive-behavioral tasks. The supragenual ACC is 
activated by pictures of spiders in phobic individuals 25 , self- 
administered but not externally- administered noxious stimu- 
lation 26 , the experience of vicarious pain 27 , painful emotions 
such as social exclusion 28 , and placebo analgesia 2 ', suggesting 
that it is involved in the cognitive and emotional modulation of 
threatening bodily stimuli. Consistent with this hypothesis, 
emotions such as anxiety and sadness intensity the negative 
hedonic quality of experimentally-induced pain, a phenom- 
enon that is correlated with hemodynamic activity of the 
supragenual ACC 30 . Further, approximately 75% of severely 
ill MDD patients treated with stereotactic lesions of the supra- 
genual ACC region (dorsal anterior cingulotomy) obtain some 
benefit from the procedure 31 . 

There are data to suggest that the stress-associated increase in 
supragenual and pregenual activity may be related to HPA axis 
function. Hubbard et al. 32 recently demonstrated that anxious 
patients with irritable bowel syndrome (IBS) show a strong 
positive coupling between BOLD response in the supragenual 
ACC and the amygdala during expectation of pain. Notably, 
administration of a CRF antagonist eliminated this positive 
coupling possibly by normalizing negative feedback by the 
amygdala on the pregenual ACC 32 . 



(2) Bilateral increases in rCBF of the amygdala were observed in 
the MDDart+ group compared with the MDDart- group (but 
not in the HCart+ versus the HCart- groups) in the ROI 
analysis (contrast A: table 2, fig. 1A). In contrast C, increased 
rCBF was observed in the right amygdala in the MDDart + 
group compared with the MDDart- group, relative to the 
HCart+ group compared with the HCart- group. An increase 
in CBF in the left amygdala in the MDDart + group compared 
with the MDDart- group relative to the corresponding differ- 
ence in healthy controls also was observed at a peak voxel 
threshold of p<0.005 (contrast C: table 3, fig. 1C). Further, 
an increase in CBF of the right amygdala in the MDDart + 
versus MDDart- group relative to controls was observed in the 
whole brain analysis with a trend towards significance for the 
left amygdala (contrast C: table 4, fig. 2C). On average, 
MDDart + patients who had higher MADRS scores showed 
significantly greater activity at the peak voxel in the right 
amygdala derived from contrasts A and C, than MDDart + 
patients with lower MADRS scores. Similarly, there was a 
trend towards a significant positive correlation between 
activity at the peak voxel in the left amygdala and 
MADRS scores in the MDDart+ group (fig. 3). The whole 
brain analysis additionally showed that rCBF was higher in 
MDDart+ patients versus MDDart- patients (but not in 
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Figure 2 Image sections from a statistical parametric map of voxel t- values which show regions where rCBF is significantly increased or decreased in 
contrasts A and C of the whole brain analysis. Stereotaxic coordinates are provided beneath each section, and are interpreted as described in the legend 
for figure 1. Abbreviations: A=anterior to anterior commissure, P=posterior to anterior commissure, L=left hemisphere, R=right hemisphere, 
spACC=supragenual ACC, pgACC=pregenual ACC, INS=insula, Tpolar C=temporopolar cortex, AMY= amygdala. Row A shows regions of the brain 
that show greater rCBF in MDDart+ versus MDDart- patients (contrast A). Al: sagittal slice showing increased rCBF to the supragenual and pregenual 
ACC (x= — 2, y=26, z=24), A2: coronal slice showing increased rCBF to the temporopolar cortex, bilaterally (x=22, y=7, z= —21), A3: sagittal slice 
showing increased rCBF to the right postcentral gyrus (x=59, y= —25, z=34). Row B shows regions of the brain that demonstrate reduced flow in the 
MDDart+ group relative to the MDDart- group (contrast A). Bl: coronal image showing decreased rCBF to the right dorsal anterior insula/ frontal 
operculum (x=36, y= — 1, z=15), B2: axial slice showing decreased rCBF to the anterior insula, bilaterally (x=36, y= — l,z=15), B3: axial image showing 
decreased rCBF to the right anterior insula (located dorsal and rostral to the cluster of voxels activated in B2 in the vicinity of the intrasulcal BA47s 
cortex 20 ) (x=26, y=20, z=6) in the MDDart+ group relative to the MDDart- group. Row C shows regions of the brain associated with greater rCBF in 
response to arterial cannulation in the MDD + group compared with the MDDart- group relative to differences between the HCart + and HCart- groups 
(contrast C). CI: axial slice showing increased rCBF to the left temporopolar cortex and right amygdala (see figure 1 legend and below for amygdala 
coordinates), C2: coronal slice showing increased rCBF to the left temporopolar cortex and right amygdala (x=24, y=3, z= — 17), C3: axial slice showing 
increased rCBF to the postcentral gyrus (x=58, y= —24, z=34). A4, B4, and C4: Key showing correspondence between image color and voxel t-value. 



HCart+ versus HCart- subjects) in mesial and rostral tem- 
poral cortex, areas that share extensive anatomical connec- 
tions with the amygdala and the medial prefrontal cortex, 
namely the right entorhinal cortex, the temporopolar cortex 
bilaterally, and the cortex along the rostral superior tem- 
poral sulcus bilaterally 33 (contrast A: table 4, fig. 2A and 
contrast C: table 6 and fig. 2C). 

Depending on the type of stressor and the specific amygdala 
nucleus examined, chronic stress causes either dendritic 
atrophy or dendritic hypertrophy of amygdala neurons in 



rodents 34 , while in humans, structural MRI is suggestive of 
decreased amygdala volumes in unmedicated patients with 
bipolar depression 35 . The functional correlates of this stress 
may be hyperexcitability of the basolateral amygdala (TSLA) 
and bed nucleus of the stria terminalis (BNST) 36 . While 
glucocorticoids usually act as a brake on hypothalamic 
CRF release, high levels of glucocorticoids potentiate the 
release of CRF from the BNST 37 and activate neurons in 
the BLA 38 . This increase in excitability of the BLA has been 
shown to occur in a timeframe relevant to this study. 
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Table 6 Whole brain analysis for contrast B showing regions where 


rCBF differed between HCart+ 


subjects 


versus HCart- subjects 


Contrast B - Whole Brain Analysis 


Increased rCBF in HCart+ vs. HCart- 


Cluster Size 


T score Z score 


p-value 


FDR p-value 


Stereotaxic Coordinates 




Region 


23 


3.67 3.58 


<0.001 


0.812 


X Y 
12 -16 


Z 
63 


R premotor cortex 


Decreased rCBF 


n HCart+ vs. HCart- 












Cluster Size 
37 


T score Z score 
3.86 3.76 


p-value 
<0.001 


FDR p-value 
0.928 


Stereotaxic Coordinates 
X Y 
-38 -41 


Z 

-8 


Region 
L parahippocampal gyrus 


Significance threshold, 


abbreviations, and interpretation of the 


coordinates are as in table 2. 









Duvarci and Pare showed that the principal neurons in the 
rat BLA increased their intrinsic activity within 1-2 hours 
of the administration of corticosterone 38 . Given that we 
have previously shown arterial cannulation to induce Corti- 
sol production and given that Cortisol increases amygdala 
excitability and CRF release, it is conceivable that the eleva- 
tion in blood flow to the amygdala in the cannulated MDD 
patients but not the HCs, is indicative of an impaired neu- 
roendocrine response to cannulation stress in the patients 
with MDD, such as an abnormal regulatory (negative feed- 
back) response to Cortisol secretion. 

Henckens et al. 35 showed that administration of hydrocor- 
tisone to healthy men acutely diminished connectivity of 
the amygdala to both limbic structures such as the hip- 
pocampus and hypothalamus as well as connectivity of 
the amygdala to cortical structures such as the middle 
frontal gyri and middle temporal gyri. The authors suggest 
that this weakening of amygdala connectivity may be an 
adaptive measure to curtail the neuroendocrine stress res- 
ponse but that continuously elevated levels of circulating 
Cortisol in MDD may result in a pattern of chronic amyg- 
dala decoupling 39 . This hypothesis receives support from a 
study that reported that healthy volunteers showed a peak 
decrease in amygdala and hippocampal activity 25-30 min- 
utes after intravenous administration of hydrocortisone 40 . 
Since [ 15 0]-H 2 0 PET imaging commenced approximately 
30 minutes after arterial cannulation (see below), healthy 
volunteers, but not depressed subjects, may have been able 
to curtail their glucocorticoid-driven stress response in the 
time period between arterial cannulation and imaging. 
Thus, a disturbance of Cortisol regulation may have con- 
tributed to the stress-associated increase in rCBF in the 
amygdalae of patients with MDD. 



(3) 



Regarding physical stressors, the amygdala shows increased 
neurophysiological response to pain in both healthy subjects 41 
and patients with MDD 42,43 . Similarly, decreased sensitivity to 
pain in PTSD is associated with a reduction in the hemody- 
namic response of the amygdala 44 . Of additional relevance is a 
recent study showing that inducing a sad mood in healthy 
subjects enhances pain sensitivity in response to a noxious 
stimulus. Greater subjective unpleasantness of the pain was 
correlated with activity in the amygdala bilaterally 45 . 
Nonetheless we note that we cannot draw any conclusions 
about the degree of physical discomfort experienced by the 
MDD versus the HC groups in this study. 
The positive correlation between rCBF in the right (and to a 
lesser extent, the left) amygdala is consistent with the finding of 
Abercrombie et al. who reported that regional cerebral metabol- 
ism in the right amygdala of depressed subjects was positively 
associated with negative affect measured using the Positive and 
Negative Affect Schedule (PANAS) 46 . Our result is also consist- 
ent with our previous finding of a significant correlation between 
Hamilton Depression Rating scores and rCBF in the left amyg- 
dala in patients with MDD 47 . Our study goes one step further by 
demonstrating that a physical stressor may induce this effect 
since the correlation between rCBF in the amygdala and depres- 
sion severity was not observed in the MDDart- group. 
A bilateral increase in rCBF was observed in the posterior 
insula (PI) in MDDart + patients compared with MDDart- 
patients (contrast A: table 2, fig. 1A; contrast C: table 3, 
fig. 1). In contrast, decreased rCBF of the left and right dorsal 
anterior insula (AI)/ frontal operculum was observed in the 
MDDart+ versus the MDDart- group (contrast A: table 2, 
fig. IB; table 4, fig. 2). 

The insula is functionally and cytoarchitecturally differen- 
tiated along rostro-caudal lines 48 . While the granular PI is 



Table 7 | Whole brain analysis for contrast C showing regions of increased rCBF in the MDD sample versus the HC sample in the high 
versus low stress condition" 

CONTRAST C - Whole Brain Analysis 



Increased rCBF in MDD vs. HC 



Cluster Size 


T score 


Z score 


p-value 


FDR p-value 


Stereotaxic Coordinates 


Region 












X Y Z 




38 


3.66 


3.57 


<0.001 


0.697 


59 -25 34 


R postcentral gyrus 


44 


3.64 


3.55 


<0.001 


0.697 


-34 5 -22 


L temporopolar cortex 


19 


3.30 


3.23 


0.001 


0.697 


24 3 -17 


R Amygdala 



Significance threshold, abbreviations, and interpretation of the coordinates are as in table 2. 
a. No region showed statistically significant decreases in rCBF for contrast C. 
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Table 8 Whole Brain Analysis showing regions of increased rCBF in the MDD sample versus the HC 
(contrast D)° 


sample in the low stress condition 


CONTRAST D - Whole Brain Analysis 


Increased rCBF in MDDart- vs. HCart- 


ClusterSize T score Z score p-value FDRp-value Stereotaxic Coordinates 


Region 


X Y Z 

14 3.41 3.34 <0.001 0.972 29 -32 4 


R Posterior Hippocampus 


Significance threshold, abbreviations, and interpretation of the coordinates are as in table 2. 
a. No regions showed statistically significant decreases in rCBF for contrast D. 



closely linked to the somatosensory cortex and parietal oper- 
culum, and provides an interoceptive map of the physiological 
state of the body, the agranular AI is thought to provide more 
abstract emotional information that flows from this corporeal 
representation. Specifically, Craig 48 postulated that sensory 
information from the PI is initially re-represented in the ipsi- 
lateral AI, but then undergoes a second-order remapping in 
the right AI, allowing for the subjective emotional experience 
associated with physiological arousal/stress. That is, the AI 
appears to play a role in translating interoceptive states into 
subjective feelings. 

Consistent with this hypothesis, the AI is activated not only by 
aversive physical stressors but also by "noxious" psychological 
stimuli such as watching a loved one in pain 49 , imagining 
oneself to be in pain 50 or believing that one is being exposed 
to harmful mobile phone radiation 51 . Greater activation of the 
right AI is also observed in response to painful stimulation 
after sad mood induction 52 . 

Since all of these studies reporting increased activity in the AI 
in response to noxious stimuli were conducted on healthy 
volunteers, it could be argued that the decrease in AI blood 
flow observed in our MDDart + sample compared with the 
MDDart- sample is reflective of an abnormal insular response 
to physiological stress in MDD. This hypothesis is consistent 
with the finding of Brody et al. 53 who reported that MDD 
patients showed increases in glucose metabolism in the left 
anterior insula (x y z coordinates: —46, 20, 8 versus a second- 
ary peak in our study at —40, 18, 6) after treatment with 
interpersonal psychotherapy. (Nevertheless, it should be 
noted that these coordinates actually localize to the frontal 
operculum associated with the inferior frontal gyrus, in the 
stereotaxic atlases of Mai et al. 54 and Talairach and Tournoux 55 
- see discussion, below). Similarly, a recent meta-analysis of 
the MDD literature reports decreased activity of the left AI 
( — 36,7,4) during exposure to negative stimuli 56 . Conceivably, 
the discrepant findings regarding functional activity of the AI 
in depressed subjects as opposed to healthy controls also may 
be related to HPA axis function. Hubbard et al. 32 found that 
compared with healthy controls, anxious IBS patients dis- 
played a strong negative correlation between the resting 
BOLD response of the AI and the amygdala during expecta- 
tion of abdominal pain. However, after administration of a 
CRF antagonist, the negative coupling between the AI and 
the amygdala was attenuated. 

In contrast, other studies of depressed patients with BD have 
demonstrated increased glucose metabolism of the right AI 57 
(34, 2, 2) increased rCBF in the right AI after induction of 
sadness 58,59 (38, 4, 0 and 44, 16, 4), and a decrease in glucose 
metabolism of the both the AI and PI after successful treat- 
ment with paroxetine 60 ( — 36, 2, 0; 40, 22, — 2) and fluoxetine 61 
(-22, 22, 6 and 36, 14, 8). Similarly, MDD patients have been 
reported to show post-treatment increases in glucose meta- 
bolism in the ventral AI (30, 18,-8 and -16,18,-16) and the 



midanterior AI (30, -2, 8 and -26, -4, 8) 62 . Consistent with 
these data, increased activation of the PI, bilaterally, and the 
left AI ( — 24, 24, 6) in response to noxious heat stimulation in 
MDD was reported 63 . Moreover Strigo et al. 43 reported that 
during anticipation of painful stimuli relative to non-painful 
stimuli depressed subjects displayed greater activation than 
controls in the right AI, although the region implicated actu- 
ally localized to the inferior frontal gyrus (44, 22, 9). 
While most studies show acutely increased activity of the AI in 
response to negative stimuli, one possible explanation for the 
decrease in rCBF in the AI in the MDDart+ versus the 
MDDart- group is that the decrease in AI activity is indicative 
of impaired sensory integration under stress. Strigo et al. 
reported that compared with healthy controls, unmedicated 
patients with MDD showed a decreased BOLD response in the 
vicinity of the right AI located in lateral orbital cortex (BA471; 
Ongur et al. 20 ) (37, 24, —7) in response to cues predicting a 
shift in the intensity of a stimulus from non-painful warmth to 
painful heat and vice versa 64 . This cortex forms part of the 
"orbital" or "sensory" prefrontal network described by 
Ongur et al. 65 which plays a major role in sensory integration. 
Strigo et al. speculated that this deficit prevents the MDD 
patients from attributing accurate emotional value to physio- 
logical safety cues 64 . Interestingly, the AI plays a role in fear 
extinction. At least 3 studies have reported increased activity in 
the AI during extinction learning in healthy individuals 66,67 . 
Gottfried and Dolan reported neural activation of the AI (39, 
12, 9 and —33, 15, 12) during the extinction-learning phase of 
an olfactory conditioning task 66 , while hemodynamic res- 
ponse of the AI ( — 35, 14, 3 and 31, 20, 5) was associated with 
extinction to conditioned shock responses 68 . In a more recent 
study, volunteers showed an increased hemodynamic res- 
ponse in the vicinity of the AI, located in the frontal oper- 
culum (-40, 20, 2 and 40, 20, -2) but a decreased 
hemodynamic response in the amygdala during the extinction 
phase of a cued fear- conditioning task 67 . Conceivably, the 
changes in rCBF in the AI reported here may be related to 
the suboptimal extinction of fear responses in patients with 
mood and anxiety disorders 69 . 

In sum, the literature on insula function in depression is con- 
fusing perhaps because the insula is a large and heterogeneous 
brain region with a highly complex pattern of internal and 
external connections. The studies that have implicated the 
AI in depression have generally reported increased activity 
in a more ventral region of the AI than reported here and this 
may be one reason for the discrepancy of our data with some of 
these aforementioned reports. The coordinates of our peak 
voxels (tables 2 and 5), are not centered in a region that is 
classically thought to be agranular AI. Rather according to the 
anatomical atlases of Mai et al. 54 and Talairach and 
Tournoux 55 , the peak voxels in our study are found in the 
dorsal insula/frontal operculum and in BA47s (Walker 12o 
in the monkey) 20 . Our group has recently shown the dorsal 
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Figure 3 Correlations between the MADRS scores and the normalized activity at the peak voxels in the left amygdala [x=— 24, y=l,z= — 20] and the 
right amygdala [24, 2, —20] in the MDDart+ group. The peak voxel for the left amygdala was obtained from the ROI analyses (contrast A, table 2, 
contrast C, table 3). Note: The normalized activity is unitless because it is a ratio of activity at the peak voxel to mean activity over the whole brain. The 
solid line represents the line of best fit for the data. 
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insula/frontal operculum to be activated specifically in res- 
ponse to interoceptive (as opposed to exteroceptive) stimuli 
at almost identical stereotaxic coordinates to those reported 
here ( — 37, — 3, 16 versus 36, — 1, 15 and 37, —3, 16 versus 
—40, —7, 17)™. Moreover Avery et al. reported that compared 
with healthy volunteers, depressed MDD subjects showed a 
reduction in BOLD activity in response to the eliciting of 
interoceptive stimuli (stomach, heart, and bladder) that was 
correlated with depression severity 70 . Clearly, the presence of 
an arterial cannula is a salient interoceptive stimulus. 
Conceivably, patients with MDD are impaired at assigning 
emotional value (including physiological safety cues) 64 to 
interoceptive stimuli, thus accounting for the decrease in 
rCBF in the dorsal AI in the MDDart+ versus the MDDart- 
groups, but not in the HCart+ versus the HCart- groups. 

(4) The MDDart+ patients displayed decreased rCBF compared 
with MDDart- patients in the sgACC (contrast A: table 2, 
fig. 2A). 

The sgACC has been implicated in pain as well as mood dis- 
orders 71 . While increased sgACC activity is usually associated 
with MDD 71 , increased activity of this region is also associated 
with pain relief. For example, activation of the sgACC has been 
reported in conjunction with placebo-induced analgesia 72 , 
deep brain stimulation-induced relief of phantom limb pain 73 
and habituation to experimentally-induced pain 74 . On the 
other hand, Lopez- Sola 63 found that compared with healthy 
controls, MDD patients failed to deactivate the sgACC and the 
pregenual ACC during painful stimulation. Conceivably, the 
cannulation-associated decrease in rCBF in the sgACC in our 
sample of individuals with MDD, is reflective of a disturbance 
in visceromotor network homeostasis which retards extinc- 
tion of somatic arousal and cannulation-associated stress. This 
hypothesis is congruent with a recent report that examined the 
neural correlates of courage as indicated by the ability to over- 
come fear in people with snake phobia 75 . Allowing a live snake 
to be brought closer to their heads without succumbing to fear 
was associated with increased hemodynamic response of the 
sgACC together with decreased activity of the amygdala. 
Moreover, the increase in sgACC activity was associated with 
a decrease in somatic arousal as measured by skin conduc- 
tance responses. 

(5) Increased rCBF in the right somatosensory cortex and precu- 
neus was observed in the whole brain analysis of the MDD 
group in the high versus low stress condition (contrast A: 
table 4, fig. 2A and contrast C: table 6, fig. 2C). 

In the meta-analysis of Apkarian et al. 18 46 out of 68 studies of (6) 
healthy subjects showed significant changes in activity of the 
somatosensory cortex during pain induction. The cluster of 
increased rCBF in the right postcentral gyrus is consistent with 
left-sided arterial cannulation. Somatosensory cortex I (SI) has 
a somatotopic arrangement and is involved in localization of 
pain on the body surface. SII, however, contains fewer noci- 
ceptive neurons and may play more of a role in spatially direc- 
ted attention towards the painful stimulus rather than sensory- 
mechanical discrimination of the stimulus 76 . SII may also play 
a role in the processing of more abstracted pain stimuli. SII is 
activated by vicarious observation of trauma 77 and hypnosis- 
induced pain 78 . In fact, a recent study even reported activation 
of SI during the processing of pain-related words 79 . Finally, the 
posterior insula, which also showed increased rCBF in 
MDDart+ patients compared with MDDart- patients, con- 
tains a major sensory association cortex 80 . 
Thus, the increased arterial cannulation-associated rCBF of the 
right postcentral gyrus in our MDD group is consistent with 3 
non-mutually exclusive possibilities: increased physical sensitiv- 
ity to the invasive stimulus in MDD, a heightened attentional 



awareness of the arterial line in our MDD sample compared 
with the HC sample, and greater physiological stressor-driven 
emotional processing in MDD. The higher rCBF in the soma- 
tosensory cortices in the postcentral gyrus in MDD subjects 
versus controls under arterial cannulation stress appears par- 
ticularly noteworthy in light of the "somatic marker hypothesis" 
of emotion proposed by Damasio and his colleagues. This hypo- 
thesis proposes that the visceral or bodily reaction that normally 
accompanies emotion or unease (the somatic marker) is experi- 
enced in a network subserved by both the somatosensory cor- 
tices and portions of the ventromedial prefrontal cortex that 
include the perigenual ACC and medial orbitofrontal cortex 81 . 
If so, then over-activity of this system in MDD in response to 
stressors may hypothetically produce the chronic sense of 
"unease" that is a frequent component of depression. 
Similarly, the precuneus has been associated with pain proces- 
sing in a number of recent imaging studies perhaps because of 
its role in self-referential processing 82 . The default mode network 
(DMN) refers to regions of the brain that are thought to mediate 
self-referential processing, are most active when the brain is in a 
resting state, and that are largely suspended when the individual 
is actively engaged in effortful cognitive tasks 83 . DMN regions 
include the dorsal and ventral mPFC, posterior cingulate/pre- 
cuneus, parietal cortex, and the superior temporal gyrus (STG) 83 . 
These regions have been demonstrated to be inappropriately 
active in patients with MDD, leading Sheline et al. 84 to hypo- 
thesize that patients with MDD have an impaired ability to shift 
their attention away from the internal milieu to meet external 
cognitive demands, and others to account for depression-related 
rumination with this model 85 . 

If correct, the DMN model not only explains the salience of 
negative emotion in MDD but raises the possibility that dispro- 
portionate internal bodily surveillance increases the intensity or 
unpleasantness of somatic signals in MDD. Consistent with this 
suggestion, chronic pain patients demonstrate abnormally ele- 
vated BOLD activity in the DMN when performing a simple 
visual attention task 86 and healthy subjects exposed to noxious 
stimuli show connectivity changes in the posterior cingulate 
cortex and the precuneus 87 . Thus the increase in rCBF in the 
ACC, inferior parietal lobe, precuneus, and STG in the 
MDDart+ sample compared with the MDDart- sample, but 
not in the HCart + sample compared with the HCart- sample, 
is consistent with the evidence for a proclivity towards higher 
activity of the DMN in MDD; a predisposition that may be 
instantiated under conditions of physiological stress. 
The MDDart- group showed increased rCBF in the posterior 
hippocampus compared with the HCart- group. A reduction 
in hippocampal volume is one of the most widely reported 
findings in MDD 88JJ0 . We have previously suggested that this 
MDD-associated reduction in hippocampal volume is more 
pronounced in the posterior hippocampus than the anterior 
hippocampus. Our group reported a reduction in posterior 
but not anterior hippocampal volume in unmedicated, remit- 
ted patients with MDD compared with healthy controls 91 , as 
well as a reduction in posterior hippocampal volume in 
patients with post traumatic stress disorder (PTSD) versus 
healthy controls 92 . Approximately half of the PTSD group in 
this study had co-morbid MDD. 

The rodent analogue of the human posterior hippocampus is 
the dorsal hippocampus. Compared with other hippocampal 
regions, the CA1 neurons in the dorsal hippocampus have 
been demonstrated to be particularly sensitive to ischemia 93 , 
possibly because of the dense concentration of NMD A recep- 
tors in this region 94 , which may confer vulnerability to glu- 
tamate-induced excitotoxicity. To our knowledge this is the 
first report of elevated rCBF specifically in the posterior hip- 
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pocampus in depressed patients versus healthy controls, and 
provides support to our previous neuromorphometric find- 
ings in depressed populations. 

In contradistinction to the widespread differences in rCBF 
that were found between the MDDart+ and MDDart- groups 
in both the ROI and whole brain analyses (tables 2 and 5), the 
HCart+ and HCart- groups did not show differences in rCBF 
in the ROI analysis, and showed only modest differences in 
rCBF in the whole brain analysis (table 6). These data poten- 
tially support the conclusion that the functional anatomical 
responses to physiological stress are exaggerated in individuals 
withMDD. 

There is increasing evidence of an overlap in the neural circuits 
that underlie physical and social pain/stress" 5 . Our results thus 
raise the possibility that individuals who are neurophysiolo- 
gically resilient to physical stressors may also be resilient to 
psychological stressors. Conceivably, future research could be 
aimed at evaluating the feasibility of using the neurophysio- 
logical response to a physical stressor as a marker of vulner- 
ability to mood disorders. 

A number of limitations deserve mention. This study was not 
designed to specifically induce pain in participants. We there- 
fore did not administer psychometric instruments to measure 
pain after arterial catheterization. Thus, although our MDD 
patients appeared more neurophysiologically sensitive to the 
stressor of arterial cannulation than HC subjects, we cannot 
draw any conclusions about differences in their psychological 
response to arterial cannulation with respect to the subjective 
experience of pain or distress. Moreover, the subjective experi- 
ence of arterial cannulation varies across subjects, such that 
some individuals report more pain and discomfort than others. 
Arterial cannulation was performed approximately 30 min- 
utes prior to the scan. Thus the initial discomfort experienced 
by subjects may have subsided prior to the commencement of 
scanning. We note, however, that glucocorticoids exert both 
rapid non-genomic, and slow genomic effects, both of which 
persist for a number of hours 96 . Henckens et al. 97 recently 
reported that administration of hydrocortisone to healthy 
individuals 75 minutes prior to fMRI scanning modulated 
amygdala activity. If the neurophysiological changes assoc- 
iated with arterial cannulation reported herein were at least 
partly related to glucocorticoid function, then MDD patients 
with and without glucocorticoid resistance would be expected 
to differ in their neurophysiological response to cannulation. 
The study would therefore have been strengthened by the 
administration of endocrine assays that would have allowed 
us to differentiate between MDD patients with and without 
glucocorticoid resistance. 

Our study was cross-sectional rather than longitudinal in nat- 
ure and therefore despite the fact that inclusion and exclusion 
criteria as well as severity of depression was indistinguishable 
across cannulation versus non-arterial cannulation studies, it 
is possible that the art+ and art- groups differed on unmeas- 
ured clinical or genetic traits that accounted for the results. 
The results reported here might not be generalizable to other 
types of stress. For instance, different results may have been 
obtained with exposure to transient, superficial pain which is 
commonly used in most experimental manipulation, is pos- 
sibly less unpleasant than deeper, visceral pain 98 , and may 
provoke less of an affective response in PET studies 99 . 
Finally, although our sample is very large for a PET imaging 
study the statistical sensitivity of our results was limited by the 
relatively small sample sizes of the non-cannulated samples. 
This reduced our power to detect differences between the 
participants who underwent arterial cannulation relative to 
those who did not, and reduced power to detect differential 



sensitivity to arterial cannulation stress between MDD sub- 
jects and healthy controls in contrast C. 
In conclusion, we raise the possibility that MDD patients are 
more sensitive - at least neurophysiologically - to invasive 
physical stressors than healthy individuals. This difference 
conceivably may account for the increased incidence of so- 
matic-related complaints in MDD. Perhaps visceromotor 
regions such as the ACC and insula, which evolved to process 
somatic pain, have been co-opted to process negative emo- 
tions. Thus on the neurobiological level, exposure to emotional 
pain and independent physical pain, may be akin to a "double- 
dose" of either emotional or physical pain. Further exploration 
of response to physiological stressors in MDD may enhance 
our understanding of the neurobiology of emotion. 

Methods 

Participants. The participants were drawn from several different PET neuroreceptor 
mapping studies, all of which were approved by the NIMH Institutional Review 
Board. All subjects provided written informed consent to participate. Detailed 
demographic and clinical data appear in table 1. 

Sixty-six individuals with MDD and 79 healthy controls participated in PET 
studies at the NIH. Blood flow images were acquired using [ 15 0]-H 2 0 before neu- 
roreceptor imaging studies that employed a range of radioligands, only some of which 
required an arterial input function for tracer kinetic analysis. Thus 51 MDD patients 
completed [ 15 0]-H 2 0 scanning with arterial cannulation (MDD art+) and 15 MDD 
patients completed [ 15 0]-H 2 0 scanning without arterial cannulation (MDD art-). 
Sixty-two HC subjects completed [ 15 0]-H 2 0 scanning with arterial cannulation 
(HC art+) and 17 HC subjects completed [ 15 0]-H 2 0 scanning without arterial 
cannulation (HC art-). 

Participants met DSM-IV criteria for MDD based on the Structured Clinical 
Interview for DSM-IV' 00 and an unstructured interview with a psychiatrist. The 
following exclusion criteria applied: left-handedness, history of bipolar disorder, 
treatment with psychiatric medication within 3 weeks of scanning (6 for fluoxetine), 
cigarette smoking, significant medical or neurological disorders, past head injury with 
loss of consciousness, serious risk of suicide, meeting DSM-IV criteria for substance 
abuse within the previous 6 months or substance dependence within the previous 5 
years, general MRI exclusion criteria, electrolyte disturbance, anemia, or positive 
drug, hepatitis or HIV screen. HC subjects had no lifetime history of a psychiatric 
disorder, and no first degree relative with a mood or anxiety disorder, as measured by 
the Family Interview for Genetic Studies (FIGS). There were no differences in 
inclusion/exclusion criteria for subjects receiving arterial cannulation versus those 
not receiving arterial cannulation. 

Twenty individuals in the MDD art+ group were medication naive and the 
remainder had been unmedicated for an average of 93 ± 120 weeks. Ten individuals in 
the MDD art- group were medication naive and the remaining subjects were 
unmedicated for an average of 28±9.0 weeks. This difference in time spent unme- 
dicated was not statistically significant (table 1). 

Eight individuals in the MDD art+ group and 2 individuals in the MDD art- group 
had a remote history of substance abuse (alcohol or cannabis). No individuals had a 
history of substance dependence. Comorbid disorders present in the MDD art+ 
group included panic disorder (PD) (n — 5), post-traumatic stress disorder (PTSD) 
(n=4), other anxiety disorders (n~3), and eating disorders (n— 2). Comorbid dis- 
orders found in the MDD art- group included PD (n~3), PTSD (n— 4), other anxiety 
disorders (n— 2), and eating disorders (n— 1). 

Depression severity was rated using the Montgomery-Asberg Depression Rating 
Scale (MADRS). Mean MADRS scores of the MDD art+ and MDD art- groups were 
25±8.1 and 23±5.4, respectively. 

Positron Emission Tomography. The MDD art+ and HC art+ participants 
underwent right antecubital venous cannulation and left radial artery cannulation 
prior to PET imaging. The latter procedure was performed by an anesthesiologist in 
the post-anesthesia care unit of the NIH Clinical Center approximately 30 minutes 
prior to commencement of the PET scan. EMLA cream was applied and a local 
anesthetic was injected into the wrist prior to arterial catheterization. The MDD art- 
and the HC art-subjects also underwent right antecubital venous cannulation 30-60 
minutes prior to the start of the PET scan. 

Each subject underwent two 40-second emission scans while resting with eyes 
closed after bolus IV injection of 10 mCi of [ 15 0]-H 2 0. 

Images were acquired in 3D mode on a GE Advance scanner and reconstructed to a 
3D spatial resolution of 6.5 mm full-width-at half-maximum. Corrections included 
measured attenuation. The [ l5 0]-H 2 0 arterial time activity curve was measured with 
an automated blood counting system which sampled blood at 3.8 ml/min during the 
[ 15 0]-H 2 0 scans in the cases who had arterial cannulae. So that the image data 
remained comparable between the subjects with versus without arterial input func- 
tions, the images of regional tissue radioactivity concentration, which are propor- 
tional to rCBF 9 , were compared after normalizing the data to exclude nonspecific 
global effects on tissue [ 15 0]-H 2 0 uptake. 
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Regions of Interest. We expected to observe MDD- associated changes in rCBF in the 
high-stress condition in the amygdala, ACC and insula. These regions were identified 
using the Wake Forest University PickAtlas. The ROI were used to constrain the 
voxel-wise analysis using the Small Volume Correction option within SPM5. 

Statistical Analysis. Analysis was carried out using SPM5. The two [ 15 0]-H 2 0 scans 

were averaged (in order to improve signal-to-noise ratio) and coregistered to a whole 

brain MR image acquired using an MP-RAGE sequence {TE— 4.94 msec; 

TR— 11.6 msec, prep time — 725 msec; delay time— 1400 msec; voxel 

size — 0.85X0. 85X 1.2 mm). After coregistration, PET images were spatially 

normalized to the standard Tl template and smoothed using an 8 mm Gaussian 

kernel. 

Analysis was performed using a general linear model: a 2X2 ANOVA with dia- 
gnosis (MDD versus HC) and arterial line (yes versus no - that is, high versus low 
stress condition) as main factors, and age as a covariate. Results were normalized by 
global activity using the SPM5 default grand-mean scaling value of 50. 

In order to measure if arterial cannulation differentially affected the rCBF of MDD 
and HC subjects, we subtracted the effects of arterial cannulation in the HC group 

(HCartH HCart-) from the effects of arterial cannulation in the MDD group 

(MDDart+ - MDDart-). In mathematical terms: [(MDDart+ - MDDart-) - 

(HCartH HCart-)]. We have labeled these contrasts as A-B=C. In addition, we 

contrasted the MDDart- vs. the HCart- group in order to evaluate group differences 
in rCBF that were specifically related to diagnosis (contrast D). 

For the ROI analyses the coordinates for peak voxel t-values exceeding a threshold 
of p<0.001 were reported if they were located within a cluster of contiguous voxels for 
which p^O.01 that was at least 100 mm 3 in size (13 voxels). To reduce Type II error 
whole brain analyses were performed post hoc to locate differences in CBF located 
outside the regions of primary interest. For the whole brain analysis, the coordinates 
for peak voxel t-values exceeding p<0.001 were reported if they were located within a 
cluster of contiguous voxels for which p<0.001 that was at least 100 mm 3 in size. 

In order to test whether the severity of depressive symptoms affected rCBF in MDD 
patients, we correlated the activity at the peak voxel of each cluster reaching statistical 
significance in contrast C (interaction) with the MADRS score for that individual. 
Specifically, the ROI analysis yielded two significant clusters with peak voxels in the 
right amygdala [24, 2, —20] and supragenual ACC [0, 27, 28]. Three clusters in the 
right postcentral gyrus [59, —25, 34], left temporopolar cortex [ — 34, 5, —22] and 
right amygdala [24, 3, — 17] reached statistical significance in contrast C in the whole 
brain analysis. However, the activity at the voxel at coordinate [24, 3,-17] in the right 
amygdala was not evaluated for correlation with the MADRS scores because of its 
close proximity to the peak voxel coordinate in the ROI analysis [24, 2, —24], which 
would have rendered the analysis redundant. We previously reported an association 
between glucose metabolism in the left amygdala and depression severity 47 , and 
therefore we also conducted a correlation analysis with the peak voxel the left 
amygdala that reached statistical significance in contrast A [ — 24, 1 , — 20] . Because of 
the a priori evidence for a positive association between amygdala activity and severity 
of depression, we conducted the correlation analyses involving the left and right 
amygdala with a one-tailed Pearson's correlation coefficient test. All other correla- 
tions were conducted using a two-tailed Pearson's test. 

The distribution of gender across the 4 subject groups was measured using the chi- 
squared test (two-tailed, p<0.05). An ANOVA was used to test for differences in age 
across the 4 groups (two-tailed, p<0.05). Independent sample t-tests (two-tailed 
p<0.05) were used to compare the MDD art + and MDD art- groups on the following 
variables: number of weeks medication free, age at illness-onset, duration of illness, 
and MADRS, scores. 
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